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The problem of dark energy is briefly reviewed in both theoretical and observational aspects. 
In the theoretical aspect, dark energy scenarios are classified into symmetry, anthropic principle, 
tuning mechanism, modified gravity, quantum cosmology, holographic principle, back-reaction and 
phenomenological types. In the observational aspect, we introduce cosmic probes, dark energy 
related projects, observational constraints on theoretical models and model independent reconstruc- 
tions. 

I. INTRODUCTION 

The Universe probed by gravitational methods looks very different from that appears optically. Thanks to a large 
number of observational windows on the gravitational side, the difference can be classified and explained by two 
additional components in our Universe: a component behaves like cold particles dubbed cold dark matter, and a 
component without significant dynamics dubbed dark energy. Both components bring big challenges to the current 
understanding of physics. 

Here we focus on dark energy [l|. Whereas observational evidences of dark energy emerges within the past 15 years, 
theoretical considerations of dark energy began almost a century ago. Those theoretical considerations are known as 
the "cosmological constant problem" . 

In 1917, the early days of general relativity, Einstein added a cosmological constant to his field equations to 
obtain a static Universe ^ . Soon in the 1920s, Pauli realized that the quantum zero point energy of a radiation field is 
too large to gravitate. Pauli's statement is already surprisingly close to the modern understanding of the cosmological 
constant problem, considering the fact that now we understand fundamental physics much better than a century ago. 

After that, the cosmological constant problem was largely ignored (or it was too hard to push forward) for a long 
time until the 1960s. In 1967, Zeldovich considered the cosmological constant problem as a fine tuning problem 0]. 
This is known as the old cosmological constant problem. Since then, a lot of early efforts have been made to eliminate 
the too large cosmological constant, as reviewed in 4]. 

In 1998, supernovae observations eventually shed light on this hard problem 01 0j and in an unexpected way. Data 
suggests that dark energy, an energy component probably coming from the cosmological constant, dominates the 
current Universe. The percentage of dark energy density is around 70% ~ not small nor too large. The Universe is 
now accelerating due to the presence of dark energy. 

In this essay we review the problem of dark energy, a.k.a. the cosmological constant problem. We first review 
the problem in the viewpoint of theoretical physics in logic order, from the aspects of symmetry, anthropic principle, 
tuning mechanisms, modified gravity, quantum cosmology, holographic principle, back-reaction and phenomenological 
models. Then we review the observational topics, including cosmic probes, dark energy related projects, observational 
constraints on theoretical models and model independent reconstructions. 



He removed this term from his equations later in 1931, due to the discovery of cosmic expansion. 
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II. THEORETICAL ASPECTS 



The theoretical challenge for dark energy can be observed from the Einstein equations 



(1) 



where A is the bare cosmological constant, as an input parameter from the theory. 

On the right hand side, the energy-momentum tensor T^i, is to be determined by a given matter configuration. 
Note that a perturbative matter state is built up from the vacuum state. Also our Universe looks not very far from 
the vacuum state. Thus to see the problem, it suffices to consider a vacuum stress tensor. 

For the vacuum, Lorentz invariance forces the stress tensor to take the form 



where Pvac = ~Pvac- Thus Pvac behaves as a renormalization of A. The quantum correction to the vacuum energy 
density can be calculated as (taking a bosonic field as an example) 



where V is the volume of space and in the last step we have used a UV cutoff A to evaluate the divergent integral. If 
an effective field theory were applicable up to the energy scale A, and natural, we could expect the following facts: 

(1) The energy density pvac gravitates. This is hard to avoid in a local field theory because the current test of 
equivalence principle has been far more precise than the quantum contribution to the mass of matter. And a local 
graviton cannot observe a difference between a vacuum loop and a loop diagram with external legs. 

(2) The effective field theory contribution to energy density pvac is not canceled by the high energy modes above 
the cutoff. This is again hard to avoid unless there are some nontrivial (and yet unknown) UV-IR correlations. 

(3) The energy density pvac is not canceled by contribution from other fields. If there were supersymmetry (SUSY), 
the bosonic vacuum energy is actually canceled by its fermionic counterpart. However SUSY, if exists, is broken on 
energy scales below lOOGeV. Thus below lOOGeV there is no known cancellation between fields. 

(4) The energy density Pvac, together with contributions from other fields, are not canceled by the bare energy 
cosmological constant A. Those two do not seem to be able to cancel each other unless gravity is modified to fit 
matter in a more dynamical way. 

Given the above "facts" , our observed vacuum energy density should be of order pvac, which is 10''^GeV'* if A = Mp, 
or lO^GeV* if A = lOOGeV. However, the dynamics of our Universe suggest a vacuum energy density pA cannot be 
much greater than the critical energy 10~^^GeV^. This is either 117 or 52 orders of magnitude smaller than the 
theoretically natural one. This is the old cosmological constant problem: 

(a) Why the cosmological constant is small? 

Moreover, recent observations show there indeed exists a component of dark energy, with energy density indeed 
of order 10~^^GeV^. The observation made the old cosmological constant problem harder to solve. Because a small 
yet non-zero number is even more difficult to explain than an unexpected zero. Also, as dark energy has different 
dynamics (dilution rate from cosmic expansion) from matter or radiation, another problem emerges: 

(b) Why dark energy density is the same order as matter density today? 

The cosmological constant problems (a) (b) are so hard that some people gave up a fundamental explanation of the 
cosmological constant and treat it as an environmental variable; On the other hand, the problem is also so inspiring 
that numerous models have been built seeking for an explanation. In the remainder of this section, we shall overview 
eight classes of those attempts. 



Symmetry has been promoted to a first principle in modern physics. It could not be better if the cosmological 
constant problem could be solved by a symmetry principle directly. SUSY Q gets near in this aspect But 
unfortunately SUSY has to be broken on large distance scales. 

Nevertheless, there are a few proposals to apply symmetry principles to the cosmological constant problem. For 
example, 't Hooft and Nobbenhuis (9| proposed that one could consider the following transformation for coordinate 
and momentum: 




(2) 




(3) 



A. Symmetry 



X' 



(4) 
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A massless scalar field or an Abelian gauge field has this symmetry. On the other hand a cosmological constant is 
forbidden. However, the mass of a scalar is not allowed either (a positive mass squire transforms into a tachyon). 
Also non-Abliean gauge theory does not have such a symmetry. There could be issues about boundary conditions in 
quantum mechanics after the above transformation as well. 

Instead o^transforming coordinate, one could as well change the sign of the metric. In a theory with g^^ — ;> "ff^i^ 
symmetry jljJl, the cosmological constant is absent. Unfortunately it is hard to realize this symmetry in four dimension 
spacetime. As another attempt [ij, one could as well scale the metric instead of changing its sign. But there is an 
instability in such a theory. 

As another class of approaches, Kaplan and Sundrum [l2j considered a counter term for the cosmological constant 

as 

C = V^J -^R - A + A„attcr(V', D^) - rn.attorlV', D I.) . (5) 

where the Lagrangian for matter field ij) and ijj have the same function form. Again there is an issue of ghost instability 
considering the and ■0 fields are coupled by gravity. 



B. Anthropic principle 

It is hopeless, if not impossible, to explain every parameter in nature from the first principles. For example, one 
cannot hope to calculate the temperature on the earth from any theory of everything. There are a vast number of 
other planets in our Universe, whose surface temperature is hugely different from that of the earth. And from the 
viewpoint of a theory of everything, the earth is not too special from any other planets to have a special temperature. 

However, a natural explanation of the earth temperature emerges, from the fact that there are intelligence livings 
on the earth, who could ask any question. This fact makes the earth special. Thus the temperature should fit for 
intelligent beings. This explanation is known as the anthropic principle [13] (of the weak sense due to Weinberg's 
classification [4]). Such an anthropic explanation of a parameter in nature makes sense if: 

(1) The parameter could vary in different parts of the Universe. On the one hand, this makes room for an anthropic 
explanation. On the other hand, this implies the value of the parameter can not be calculated from the first principles 
(although its distribution may be calculated). 

(2) The existence of intelligence is sensitive to this parameter, and the observed value (at the time of observation) 
prefers such an existence. 

(3) Additionally, to relate an intelligent observer to us humans, one has to assume we are typical observers, and 
propose a measure to define the typicality in a probably infinite Universe. 

The above are the general considerations of the anthropic principle. Those considerations may be applied to the 
cosmological constant problem. This is because: 

(1) There might be different cosmological constants in different parts of the Universe, in space scopes beyond 
our observable Universe. This statement is (so far) not testable by the (current) definition of observable Universe. 
However, fundamental theories suggest such a possibility. Es pec i ally , in string theory, in the past decade a landscape 
of (semi-classical) meta-stable vacua solutions are discovered [la] [1^] . The number of such vacua in string theory may 
be of order 10^°° [13], considering one could use different Calabi-Yau manifolds to compactify the extra dimensions, 
and add different branes and fiuxes [l^. Moreover, the theory of eternal inflation [l^, as well as the quantum 
mechanical many world conjecture, open up possibilities to populate the string landscape. If the above picture holds, 
the cosmological constant indeed ends up as an environmental variable. 

(2) The existence of intelligence is indeed sensitive to the cosmological constant. This is because a large and positive 
cosmological constant will dilute the structure of the Universe before useful structures (like galaxies) are formed 2^ ■ 
Whereas a large (in absolute value) and negative cosmological constant will drive an early collapse of the Universe. 
However, it is still too early to conclude that our observed value at the time of observation is most suitable for 
intelligence. For example, there is a Boltzmann brain problem ^21j, indicating that a (nearly) pure de Sitter Universe 
may be more suitable for intelligence, of a thermal fluctuation type. 

(3) Finally, it is in debate if we can typically consider ourselves as typical observers, and consider the observed 
Universe as a typical one in the multiverse. Moreover, the measure problem (2^ is still severe. Currently several 
measures coincide with each other and a few other measures are ruled out. Thus it seems that different efforts converge 
in the right direction. However, no guiding principle is available so far for the measure construction thus the problem 
should not be marked as solved. 
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C. Tuning mechanisms 

Another direction for dark energy is to introduce a field, wliicli could dynamically cancel the cosmological constant, 
no matter how large it was originally. Such kind of a scenario is called tuning mechanism. For example, for the old 
cosmological constant problem, one may have expected to have a scalar field 0, which is related to the stress tensor 




V2(/) (xTi; (xR . (6) 

If also depends on and drives 4> to some stable value (j)Q, the equilibrium state sd, (j) — (po will have a vanishing 
Ricci scalar thus a vanishing cosmological constant. 

Unfortunately, Weinberg j4j showed that it is difficult to realize such a mechanism. The above conditions are proven 
to be equivalent to the statement that the matter Lagrangian takes the form 

£ = e^^V^Co. (7) 

Thus unless one fine tunes Cq, there will not be a stationary point for 0. Also, if asymptotically rolls to — (X), the 
effective Newton's constant Gn vanishes, because (j>Q can be viewed as a redefinition of Gjv. Neither of the above 
possibilities is satisfactory. 

There are some recent attempts trying to get around the above problem. For example, one can consider a non-trivial 
extra dimension and then Weinberg's proof no longer applies {23 ]. However, attempts in this direction also encounter 
problems such as a singularity in the bulk, or violation of the Lorentz symmetry. 



D. Modified gravity 

Dark energy (as well as dark matter) is identified via gravitational probes. Thus once one meets difficulties in the 
explanation of dark energy, it is a natural idea to step back and question if the current gravity theory can be applied 
to cosmological scales (or even galactic scales). 

There have been a large number of theories of modified gravity, and a large proportion of which has been used 



as models of dark energy. An incomplete list of which include .f{R) gravity [251 |l26ll[27|, MOND [281 and TeVes [29|, 
DGP 30], scalar-tensor theories [31] 32], Gauss-Bonnet ^33] [34] and Lovelock [35] gravities, Horava-Lifshitz gravity 
[36], /(T) gravity jsvl ]. conformal gravity [s^, fat graviton and so on. Those scenarios open up a huge variety of 
possibilities for dark energy. 

The attempts of modifying gravity starts as early as the early 1960s, known as the Brans-Dicke gravity [3l|. In 
Brans-Dicke gravity, a scalar field is added into the gravitational action 

S^l d^x^g {^-d^R ~ ^{d^f^ . (8) 

Instead of inserting a scalar field, one could also modify the Einstein-Hilbert action by using an alternative curvature 
invariant. As an example, f{R) gravity modifies the Einstein-Hilbert action as 

" d^x^ i? ^ — ^ / d^x^/^ f{R) . (9) 



2 

To make specific predictions one has to specify /. For instance, one could choose f{R) — R — A'^/i?, where fi is a 
constant (26j . 

Along this direction, one could also construct more general curvature invariants. For example, the Gauss-Bonnet 
term 

R^B = R'"''"'R^.pa - 4i?^'^i?^, + R^ . (10) 

In 4 spacetime dimensions Rq^ is a total derivative thus does not modify the local dynamics of gravity. However, one 
could couple -Rqb to a scalar field as /(0)i?GB ^ model of dark energy [34]. 

The DGP scenario is another popular possibility of modifying gravity _3Q;i. In the DGP scenario, two gravity 
theories are considered at the same time: one located on a 3-1-1 dimensional brane, and another on a 4-1-1 dimensional 
brane, with their own Einstein-Hilbert actions (while the two gravity theories share the same metric: the 3-1-1 metric 
is the 4-fl metric restricted on the 3-brane). The existence of the 4-fl brane leads to large scale modifications of 
gravity. 
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E. Quantum cosmology 

The cosmological constant problem is a problem from the energy of quantum fluctuations, which gravitate. Thus 
it is not surprising if the solution could have something to do with quantum gravity, even though the problem mainly 
affects extremely large scales. 

For example. Hawking [s^ji^] provided a mechanism to show that the cosmological constant is probably zero. He 
constructed a 4- form field strength, which could vary and change the cosmological constant. Then the wave function 
of the Universe indicates a probability distribution of the effective cosmological constant A (before normalization of 
the probability distribution) to be 

P(A) = e^. (11) 

Thus A — )■ from above is the most probable value of the cosmological constant when the Universe is created via 
the Hartle-Hawking wave function. Unfortunately, this scenario produces an empty Universe, with only a vanishing 
cosmological constant. On the one hand, there is no radiation or matter in this approach, and on the other hand, our 
present accelerating Universe (if not coming from another source) falsifies the mechanism. 

Later, the probability is corrected by considering the decoherence effects [4l|. In this case the Universe could be 
born at the string scale. However, in this case this mechanism is no longer very related to the cosmological constant 
problem. Instead, the mechanism is used as a program to select the original Universe. 

As another quantum approach to dark energy, Tye j43 | considered a possibility that the vacuum of our Universe 
may experience fast tunnelings through the string landscape, until settling down at a place with a small cosmological 
constant. The idea is supported by the resonant tunneling effect in quantum mechanics. Resonant tunneling is a well 
established effect in quantum mechanics. Recent development also shows evidence that resonant tunneling can exist 
in quantum field theory 0] H gi] . 

One could also relate quantum features of cosmology to dark energy from a cosmological seesaw mechanism [47l . |48| . 
It is interesting to notice that if SUSY is broken at the TeV scale, there is a hierarchy of scales 

Ma ^ Mi^sy/Mp, (12) 

And the SUSY scale could also be replaced by any scale that stabilizes Higgs mass. It is argued that the hierarchy 
could be resolved by a seesaw mechanism relating two interacting Universes. The coupling is introduced through the 
wave function of the Universe (49j . It is so far not clear yet how to introduce such a coupling at a more fundamental 
level. 



F. Holographic principle 

Effective field theory is one of the most important concepts in high energy physics. The cosmological constant 
problem, as a problem of renormalizing the cosmological constant, can also be thought of as a problem in the sense 
of effective field theory. However, the application regime for effective field theory seems not as well understood as 
people thought years ago: Not only an increase of energy scale, but also an increase of the number of relevant degrees 
of freedom, seems to lead to a break down of effective field theory. 

There are quite a few examples for the above observation. Here we only consider the most relevant one for our 
purpose: Consider the vacuum energy within a volume with radius L [5(i|. The volume under consideration should 
not form a black hole thus the effective field theory should have a cutoff Auv satisfying 

L^Afjy ^E< LMl , (13) 

where the very right hand side is the energy of a black hole of the same size. This inequality can be recast in terms 
of energy density as 

p = 3c^M^L-^ , (14) 

where c is a parameter introduced in [5l| . which is naturally a constant of order one. 

As pointed out in (5l| . when L is chosen the event horizon of our Universe, p behaves as a component of dark 
energy, dubbed the holographic dark energy (HDE) . The equation of state for this component is 



1 2Vfihdc ,^r\ 

Whdc = 5 , 15) 

3 3c 
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where dhdo is the relative energy density of HDE. Thus when HDE subdominates, it behaves as curvature; and when 
it dominates the Universe, the Universe behaves as quasi-de Sitter for reasonable c. 

The equation ()14|) can be also understood in several other ways. For example, as shown in [s^ . the Casimir energy 
in de Sitter space has a quadratic divergence. Thus the Casimir energy takes the form (|14p and de Sitter expansion 
may be self-driven. Alternatively, if one considers the event horizon of our Universe as a holographic screen, the 
acceleration of Universe behaves as an entropic force 53], or a consequence of quantum information M- Also, if 
non-perturbativc gravitational excitations behaves as a gas component with holographic dispersion relation, one could 
also arrive at (ITO [55]. 

The event horizon is a natural choice, but not the unique logical possibility. It is also proposed that the Ricci scalar 
[56j . or the time of the Universe [13], could be used as the IR cutoff as well. Also the parameter c could be a function 
of time or redshift [l!] . 



G. Back-reaction 



The Einstein equations are non-linear. The non-linearity brings a great amount of difficulties as well as challenges 
to applications of the theory of gravity. 

One of the distinguishing features of a non-linear theory is that the time evolution and the operation of averaging 
may not commute. In the conventional approach of cosmology, we first approximate our Universe to be homogeneous 
and isotropic. This is true only on large scales. We are averaging over smaller structures before considering time 
evolution of the scale factor. Thus there will be a correction from the non-commutivity of averaging and time evolution. 
The correction is dubbed back-reaction [59], and the question is if the correction could become large and accelerate 
our Universe. 

The back-reaction term appears in the cosmological equations as [60j 

3- = -4ttG{p) + Q, (16) 
a 

3H' = 87rG{p)-^{R3)~^Q, (17) 



dt{p}+3H{p}=0, (18) 

where is the spatial curvature, and Q is defined as a combination of the variances of the expansion scalar 9 and 
the shear a as 

Q^^((0-(e)')-2(0. (19) 

It is still an open question if Q could be large enough to drive the cosmic expansion. Also, in this approach, one still 
has to solve the old cosmological constant problem, to have a not-too-large original cosmological constant. 

One could also consider whether super- Hubble inhomogeneities could back-react the evolution of our Universe. This 
is a more ambitious approach. Because if there is a screening effect of cosmolo gica l constant form the back-reaction, 
one can solve the old cosmological constant problem together. One can prove [6l|[6^ that if there is only curvature 
fluctuation, the super-Hubble inhomogeneities cannot back-react onto our Hubble scale dynamics. When there is 
isocurvaturc fluctuation, the effect of super-Hubble back- reaction is still in debate (63j . 



H. Phenomenological models 



As another huge class of attempts, dark energy could as well be a simple fluid component, for example, a scalar 
field. The simplest example of this class is a quintessence field with a conventional scalar [g^, IIH] 



5 = y d^x^—g 



1 



At the homogeneous and isotropic background level, the energy density and the pressure takes the form 



(20) 



(21) 
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When the potential V is flatter than 

an accelerating solution for our Universe is admitted. In the limit when the potential V is completely flat, the scenario 
returns to a massless scalar field, plus a decoupled non-dynamical cosmological constant. A quintessence field always 
has equation of state w > —I. 

The quintessence scenario could be modified in various ways. For example, one could consider a scenario of phantom, 
which introduces a quintessence field with a wrong-sign kinetic term (66j . In this case, one get a field with w < — 1 as 
dark energy. The phantom field is unstable when coupled to matter (even gravitationally) . Even worse, in a Lorentz 
invariant theory, the decay rate of phantom is divergent [67| . Recently, it is argued that a dynamical breaking of 
Lorentz invariance due to the decay of phantom itself may avoid this problem [68| . The combination of quintessence 
and phantom fields also makes a scenario of quintom [6^. 

As another direction, one could generalize the standard kinetic term as [T^I 



^^g^-d^ipd^ip ^ P (-^g^-d.ipd,^, ifi ) , (23) 



where P is a general function of its two arguments. This generalization is known as k-essence. The scenario can be 
further generalized to ghost condensation [71,, J2^ and kinetic gravity braiding [t^ . 



The field theoretic approach of dark energy can also be generalized to spinor fields [TJ] , vector fields [75[ , and p-forni 
fields [t^. Also, dark energy could be particles with exotic properties [55|[77j. 

As another class of phenomenological models, dark energy may as well be a fluid with a modified equation of state. 
For example, Chaplygin gas [78|, (73] is a model in this direction. Chaplygin gas, originally introduced to model the 
air fiow around the wing of an aircraft, is a fluid with equation of state 

p = -A/p . (24) 

Later, the scenario is generalized to [13] 

P = -A/p" (25) 

to better fit the data and to unify dark matter and dark energy. 

As another example of fluid dark energy, dark energy might be a fluid with viscosity f sij 

Ti_,„ = pUf^u^ + {p- 3HC){gi_,^ + u^Uy) . (26) 

There are two generic features for most phenomenological models of dark energy. Theoretically, one still need to 
worry about the cosmological constant problem because the large vacuum zero point energy is not canceled in a 
natural way. Phenomenologically, one should consider dark energy perturbations as well in those models, because 
dark energy has local degrees of freedoms in those cases. The latter could be a window to distinguish field, fluid, 
particle scenarios of dark energy between a simple cosmological constant. 



III. OBSERVATIONAL ASPECTS 

A. Cosmic probes of Dark Energy 

The most common approach to probe dark energy is through its effect on the expansion history of the Universe. 
This effect can be detected via the luminosity distance ^^(z) and the angular diameter distance Da{z). In addition, 
the growth of large-scale structure can also provide useful constraints on dark energy. In this section, we briefly review 
some mainstream cosmological probes of dark energy. For a more detail review on current mainstream observational 
probes, see Ref. [82j . 



1. Type la supernoave 



A white dwarf star can accrete mass from its companion star; as it approaches the Chandrasekhar mass, the 
thermonuclear explosion will occur, representing itself as a type la supernoave (SNIa). This mechanism implies that 
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SNIa can be used as standard candles to measure the luminosity distance dL^z), which takes the form ^^(z) = /p -jjfp^- 
Clearly, (Il^z) contains the information of the cosmic expansion history H{z) through an integration. 

In 1998, using the SNIa data from the Hubble space telescope (HST) observations, Riess et al. 5] first discovered 
that the expansion rate of our Universe is accelerating. Soon after, the result was confirmed by Perlmutter et al. Q 
based on the analysis of 42 SNIa. The discovery of the Universe accelerating expansion was named "Breakthrough 
of the Year" by Science Magazine in 1998. After that, the surveys of SNIa has drawn lots of attention. Famous 
teams focusing on this field include the High- 2; team ^8^], the Supernova Legacy Survey (SNLS) fs^, the ESSENCE 
[ssj l. the Sloan Digital Sk y S urvey (SDSS) Supernova (SN) Survey iSfii], and so on. Recent SNIa dataset includes the 
three-year SNLS dataset [83] consisting of 472 SN and the Supernovae Cosmology Project (SCP)'s Union 2.1 dataset 
[ssj consisting of 580 SN. When combined with other cosmological probes, these data could reduce the (Ict) error of 
the dark energy EoS to less than 0.1. 

For a SNIa, once its redshift Zi, absolute magnitude Mi and apparent magnitude rrii are determined, one can 
construct the theoretical (assuming specific dark energy) and observational distance modulus, 

lJLobs,i ^ mi - Mi, ijLth,i^^^ogiodL{zi) + 2b, (27) 

and the for the SNIa data can be constructed as Xsat — X^i ■ In this procedure, people often analytically 

marginalize the nuisance Hubble constant Hq (89| . 

Due to the absence of the knowledge of its detailed explosion mechanism, SNIas are not intrinsically standard 
candles. Other systematic errors of the SNIa come from the photometry, the identification of SNIa, the host-galaxy 
extinction, the gravitational lensing, and so on. These systematic errors are now the major factor that confines the 
precision of SNIa data. To enhance the precision, improvements on the photometric technique, as well as better 
understandings of the dust absorption and the SN explosions, are needed. 



2. Cosmic Microwave Background 

As the legacy of the cosmic recombination epoch, the cosmic microwave background (CMB) contains abundant 
information of the early Universe. Current CMB experiments include the well-known Wilkinson Microwave Anisotropy 
Probe (WMAP) the Planck satellite fgH], and so on. 

From the CMB observations, some distance priors can be extracted to constrain dark energy. For example, the "shift 
parameter" [92| defined at the decoupling epoch z* takes the form i?(z*) = r2miJo(l + z^)Da{z^). This parameter 
nicely describes the full CMB leverage on the cosmic expansion history, and is a good approximation for models not 
too much deviated from ACDM. Another distance ratio is the "acoustic scale" I a, defined as (1 + (''« 
is the sound horizon). This quantity represents the CMB multipole of the acoustic peak location. As high redshift 
distance indicators, R and Ia are especially useful in breaking the degeneracies between cosmological parameters when 
combined with other cosmological probes. 

The CMB data could also be used to probe dark energy through the Integrated Sachs Wolfe (ISW) effect [o^ . This 
large scale anisotropy effect is caused by the variation of the gravitational potential during the epoch of the cosmic 
acceleration, and has been detected through the cross relation between CMB and the LSS at about 4(7 (9^. 



3. Baryon Acoustic Oscillations 

Baryon acoustic oscillations (BAO) refers to overdensities or clusterings of baryonic matter at certain length scales 
(~150 Mpc in today's Universe) due to acoustic waves which propagated in the early Universe [95, 96]. It provides a 
"standard ruler" for cosmological observations, and can be measured at low redshifts z < 1 through galaxy surveys. 
Famous BAO measurements include the Two-degree-Field Galaxy Redshift Survey (2dFGRS) the SDSS j98|], the 
WiggleZ dark energy survey [9^, and so on. When comb ined with the SNIa and CMB data, the latest BAO data 
from the WiggleZ survey yields a resuh w = -1.03 ± 0.08 [Tooj . 

From galaxy observations, the BAO scales in both transverse and line-of-sight directions are obtained, corresponding 
to the quantities r{z)/rs{zd) {r{z) is the comoving distance; Zd is the "drag epoch") and rs{zd)/ H{z). Widely used 
BAO distance measurements include the A parameter, the distance Dv{z), the distance ratio rs(zd)/ Dy{z), and so 
on. Current measurements distribute at the redshift region z ~ 0.1 — 0.75. 

Since the BAO measurement only requires the determination of the three dimensional positions of galaxies, they 
are less affected by astronomical uncertainties. 
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4- Weak lensing 

Weak lensing (WL) is the slight distortions of distant objects' images, due to the gravitational bending of light 
by structures in the Universe. The mass and positions of the lens depend on the matter distribution on the light 
cone, while the distances to the objects and lens are determined by the geometry of the spacetime. So WL probes 
dark energy through its effects on both the cosmic expansion a nd g rowth history of the strucuture. The current WL 
project is the Canada- France-Hawaii legacy survey (CFHTLS) [101]. 

The effect of WL on the distance objects represent on the distortions of the shapes of the galaxies [l02j | ^, with 
typical scale ^0.01. This signal is much smaller than the typical scale of the intrinsic deviation in the galaxy shape 
(about 0.3-0.4) , and can be only detected from a large number of galaxies. The most commonly used method is 
the cosmic shear statistics. In 2005, based on the cosmic shear data of CFHTLS, Hoekstra et al. [l04| reported a 
constraint oi w < —0.8 at the 68% CL. 

Although current WL data are not powerful enough compared with SNIa and BAG, WL is a very promising probe 
to detect dark energy. The underlying WL data are rich, while the observational and statistical techniques of the WL 
probe are under continually improvement. 

5. Other dark energy probes 

Galaxy clusters: Their distribution marks the LSS. Since the comoving volume of t he Un iverse is affected by the 
geometry, the cluster abundance also encodes the information of the expansion history [l05| . 

Gamma-ray bursts: As the most luminous electr oma gnetic events in the Universe, gamma-ray bursts (GRBs) may 
help us to extend the Hubble diagram to z ~ 8 |106{ . However, due to the large errors and calibration problem, 
current GRB data are not able to provide extensive and reliable probe to dark energy. 

X-ray observations: Like the SNIa and BAO data, the X-ray fgas (defined as Mgas/Mtot) measurements can also 
probe the redshift-distance relation, with the dependence fgas oc dL{z)DA{zY'-^. In 2008, based on the Chandra X-ray 
data, Allen et al. [lo3l reported Q.k = 0.86 ± 0.21. 

Hubble parameter measurements: Precisi on m easurements of Hq can help to break the degeneracies, and are thus 
widely used in dark energy constraints. In |l08l |. Riess et al. reported Hq = 73.8 ± 2.4km/s/Mpc. Besides, the direct 
measurement of H{z) from the differential ages of passively evolving galaxies can also provide valuable constraints 

Cosmic age test: Cosmic age problem is a "smoking-gun" of evidence for the existence of dark energy ^. However, 
the existence of some extremely old globular clusters an d quasars shows that the cosmic age problem has not been 
completely removed by the introduction of dark energy |l09l . IllOj ] . In addition, the "look back time-redshift" |lll| 
data obtained from the age of old passive galaxies can be included into the analysis and constrain dark energy. 

Growth factor data: The growth factor defined by / =_d\n6 / dlna can be measured from the redshift distortions 
and Lyman-a Forest data. Current growth factor data jll2l | have large errors and are still less powerful. 

B. Dark energy projects 

According to the DETF report [ll3l |. the dark energy projects can be classified into four stages: completed projects 
are Stage I; on-going projects, either taking data or soon to be taking data, are Stage II; intermediate-scale, near- 
future projects belong to Stage HI; larger-scale, longer-term future projects belong to Stage IV. In this paper, we 
will i ntro duce some mainstream Stage IV projects. For a comprehensive list of the dark energy experiments, see Ref. 



^ It should be noted that the lensing of CMB can also be used to probe dark energy. In 2011, using the Atacama Cosmology Telescope 
lensing measurements, Sherwin et al. llOSH reported that the CMB measurement alone favors cosmologies with w = —1 dark energy 
over models without dark energy at a 3.2f7 level. 

^ In the matter-dominated SCDM model, the present age of the Universe is only to = 3^^' = 9Gyr, and contradicts with the detections 
of many older than lOGyr objects in astronomical observations. 
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1. Large Synoptic Survey Telescope 

A most ambitious ground-based dark energy survey project is the Large Synoptic Survey Telescope LSST [ll6| . 
which is an USA NSF/DOE collaboration program. Its design is an 8.4-meter ground-based optical telescope to be 
sited in Cerro Pachon of Chile. Over a IQ-year lifetime, LSST will obtain a database including 10 billion galaxies and 
a similar number of stars, studying dark energy through a combination of the BAO, SN and WL techniques. Since 
its compelling scientific capacity and relatively low technical risk, LSST was selected as the top priority large-scale 
ground-based project for the next decade of astronomy in the Astro2010 report [ll7j |. The project is scheduled to 
have first light in 2016 and the beginning of survey operations in 2018. 

2. Square Kilometer Array 

Another ambitious large-scale ground-based telescope is the Square Kilometer Array (SKA) [ll8j |. an international 
collaboration program. SKA has a collecting area of order one square kilometer and is capable to operate at a wide 
frequency range (60 MHz - 35 GHz). It will be built in the southern hemisphere (either in Australia or in South 
Africa), and the specific site will be determined in 2012. SKA will probe dark energy by BAO and WL techniques 
via the measurements of the Hydrogen line (HI) 21-cm emission in normal galaxies at high redshift. Its construction 
is scheduled to begin in 2016 for initial observations by 2020 and full operation by 2024. 

3. Wide Field Infrared Survey Telescope 

A most exciting space-based dark energy project is the Wide Field Infrared Survey Telescope (WFIRST) [ll9j | , which 
is an USA NASA/DOE collaboration program. WFIRST is a 1.5-meter wide-field near-infrared space telescope, orders 
of magnitude more sensitive than any previous project. It will determine the nature of dark energy by measuring the 
expansion history and the growth rate of large scale structure. A combination of the BAO, SN and WL techniques 
will be used. Since its compelling scientific capacity and relatively low technical risk, WFIRST was selected as the 
top priority large-scale ground-based project for the next decade of astronomy in the Astro2010 report (llTj . The 
project is scheduled to launch in 2020 and has a 10-years lifetime. 

4- Euclid 

Another exciting space-based dark energy project is the Euclid [l20j |. which is an ESA project. Euclid is 1.2 m 
Korsch telescope, with optical and near-infrared (NIR) observational branch. It will also make use of several secondary 
cosmological probes such as the ISW, CL and redshift space distortions to provide additional measurements of the 
cosmic geometry and structure growth. After 5 years' survey, Euclid will measure the DE EoS parameters wq and Wa 
to a precision of 2% and 10%, respectively, and will measure the growth factor exponent 7 with a precision of 2%. 

5. International X-ray Observatory 

The last Stage IV project is and the International X-ray Observatory (IXO) [l2l|, which is a partnership among 
the NASA, ESA, and the Japan Aerospace Exploration Agency (JAXA). It is a powerful X-ray space telescope that 
features a single large X-ray mirror assembly and an extendable optical bench with a focal length of ^ 20 m. With 
more than an order of magnitude improvement in capabilities, IXO will study dark energy thro ugh t he WL technique. 
It was selected as the fourth-priority large-scale ground-based project in the Astro2010 report jll7l |. and is scheduled 
to launch in 2021. 



C. Observational Constraints on Theoretical Models 

In this section we briefly review some research works concerning the cosmological constraints on various theoretical 
models. 
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1. Scalar field models 

As the most popular phcnomcnological dark energy models, scalar filed models can mimic a dark energy component 
with arbitrary properties (e.g., arbitrary w). So the issue of the observational constraints on the scalar field models 
is somewhat similar to the observational probe of the dynamical behavior of dark energy. 

Different classes of scalar field dark energy models have different properties. For quintessence or phantom, w > —1 
or w < — 1, while for quintom w can cross — 1 ^. Current observations still can not determine whether dark energy is 
quintessence, phantom or quintom. For example, the Seven-year WMAP measurements [l23| gave w = —1.10 ±0.14, 
with all the possibilities exist. 

Another popular issue ofte n dis cussed in the literature of the scalar field models is possible interactions between 
dark energy and dark matter |l24| . Phenomenologically, interaction can be introduced in the form 

Pm + 3Hpm = Q, Pde + 3H{pde + Pde) = -Q, (28) 

where Q is the interaction term. Current data put tight constraints on interactions between the dark sectors. In 
p~25i] ■ based on the analysis of the SNIa+BAO+CMB data, Het et al. considered the form Q — 3C,Hp, and obtained 
C ^ 0.01 for the cases oi p = p,„ and p = p„i + pi^^. Similar results were obtained in e.g. 126j. 



2. Chaplygin gas models 

The standard Chaplygin gas model with jo = A/p has been ruled out by the SNIa, BAO and CMB observations at 
more than 99% confidence level (CL) (l27| [ l28 | . 

The generalized Chaplygin gas model with p = A/ p" is once believed as a promising candidate to unify dark matter 
and dark energy (hereafter noted as "UDM" ) . The evolution of its energy density takes the form, 

p{t) ^ po[{l ~ n,) + n,a-'"^^+''^]T^ , (29) 

where Q^, is the effective matter ratio determined by A and and integration constant. The model can nicely fit the 
observations of the cosmic expansion history when a ~ 0. Especially, the case of the ACDM model corresponds to 
a = 0. 

However, it is later proved that |l29| this UDM scenario could produce oscillations or exponential blowup of the 
matter pow er sp ectrum, which is inconsistent with observations. In this way, the parameter space is confined to 
^ 10~^ jl30l |. extremely close to the ACDM case. So the generalized Chaplygin gas model as a UDM has been 
effectively ruled out. 



3. Holographic dark energy models 

In [iH, Li proposed the HDE model with the future event horizon as the IR cutoff. Hereafter, a lot of works have been 
performed testing this models using the observational data of SNIa, BAO , CM B, WL, and so on. These works all show 



that the HDE model can provide a good fit to the data. For example, in [l31j . by using the Consitution+BAO+CMB 
data, Li e t al. found that Xacdm = 

467.775 and 

Xhde — 465.912, implying that HDE is a competitive dark energy 

candidate [l32| . 

Numerical studies of the HDE model show that the data favor c < 1 [T33l - [l36j . Recent studies have obtained the 
constraints c < 0.9 at the 95% CL. The data also put tight constraints on possible interactions between HDE and 
dark matter. Generally speaking, for an interaction term with the form Q = iC^Hp (where p can be p™, Pde, Pm + Pde, 
and so on), a constraint of |C| ^ 0.01 can be obtained based on the anal ysis of the SNIa+BAO+CMB data [l37| |. 

As mentioned above, the Ricci scalar [5^ and the time of the Universe [57l | have also been proposed as the IR cutoff. 
However, later researches showed that these two mode ls are not favor ed by the data. Compared with the original HDE 
model, the xLin of these models can be 15-30 larger [ml [iM Eli] based on the analysis of the SNIa+BAO+CMB 
data. 



The value of dark energy EoS w is crucial to the future fate of our Universe. The Universe will end up writh a "big rip" if «) < — 1 [122| |. 
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4. DGP model 



Although it is a simple mechanism for the cosmic acceleration without introducing the dark energy component, 
the DGP model can not provide a satisfactory fit to the data. Firstly of all, the exp ansio n history predicted by the 
DGP model leads to appa rent inconsistency between the SNIa, BAG and GMB data |128{ . implying that this model 
is not acceptable. In (l38j . by using the Union+BAG+GMB data Rubin et al. showed that the DGP model leads to a 
^Xmin ~ 1^ compared with the ACDM mode l. Secondly, the history of the structure formation in the DGP model is 
also disfavored by observations jl39l |. In [l40l |. Fang et al. showd that this model over predicts the low-/ anisotropics 
of the GMB temperature power spectrum. Gompared with the AGDM, DGP is excluded at 4.9ct (with curvature) 
and 5.8ct (without curvature) levels. Thus, the situation of the DGP model is far from optimistic when confronted 
with cosmological observations [141i] . 



5. Other modified gravity models 

f{R) models can mimic indistinguishable expans ion h istori es as the AGDM model, so to test them we need to lend 
help from the cosmic growth history observations |l42l |. In [l43| [. Song et al. parametrized the f{R) model by the 
"Gompton wavelength parameter" B ^, and found that the model predicts a lower large-scale GMB anisotropics by 
reducing the ISW effect. So far, tightest constraints on the f{R) models come from the observations of the galaxy 
clusters. In |144| . Lombriser et al. revisited the model studied in and reported i?o < 1-1 x lO^'^ at the 2a GL, mainly 
due to the inclusion of the cluster abundance data. 

Gauss-Bonnet gravity predicts strong growth of perturbations on smaller scale s, in compatible with the observed 
galaxy spectrum, unless the deviation from the Einstein gravity is very small |l45j . Thus this model has been 
effectively ruled out. 

Based on the constraints on the PPN parameter 7 from the solar system and binary pulsar observations [l46j | [l47| , 
the Brans-Dicke theory has been tightly constrained to the level that its cosmological effects are insignificant. 



6. LTB models 



LTB models are the most commonly considered inhomogeneous scenarios to explain the cosmic acceleration without 
introducing dark energy. In the simplest class of such models, we live close to the center of a spherically symmetric 
Gpc scale void, and our local region has a larger Hubble parameter than the outer region due to the spatial gradients 
of the metric. These mo dels can successfully explain the dimming of the distant supernovae, providing a nice fit to 
the SNIa data [Hilflii], comparable with the AGDM model. 

In [150], Alnes et al. pointed out that, if it happened that we do not live precisely near (less than ~20 Mpc |151l] ) 
the center of the void, the observed GMB dipole would become much larger than that allowed by observations. So 
there is a fine-tuning problem in the void models. What is worse is that, even if we happen to live very close to the 
void center, the existence of the off-center galaxy clusters will le ad t o an observationally kinetic Sunyaev-Zeldovich 
effect due to their relative motion between the GMB frame. In jl52j |. Bellido et al. demonstrated that the limited 
observations of only 9 clusters with large error bars already rule out LTB models with void sizes greater than 1.5 Gpc 
and a significant underdensity. 

There have been many other intere sting methods proposed to test the LTB models, including the Hubble constant 
measurement |108{ . the r edshift drift 153 1, the ionized Universe mirror 154 1. the cosmic age test [155i] . the constant 
curvature condition (l56j |. and so on. The test of the LTB models is related to the test of the Gopernican principle, 
and this issue has drawn a lot of interests in recent years. 



7. Comparison of dark energy models 



The analysis alone cannot provide an effective way to make a compa rison among dark energy models. To enforce 
a comparison, a general way is to employ the information criteria (IG) [157| . The most frequently used IG are the 



5 The definition is B = 4^ ^tf/^1'"° . B > is the stable branch. 

/ dH/d In a 
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Bayesian information criterion (BIG) |l58j | and the Akaike information criterion (AIC) |l59l |. defined as, 

mC = -2\n£,nax + k\nN, AIC = -2\nC,nax + ^k, (30) 

where Lmax is the maximum likehhood (satisfying —2\nCmax = Xmin under the Gaussian assumption), k is the 
number of the parameters of the considered model, and N is the number of data points used in the fit. Glearly, these 
statistics favor models that give a good fit with fewer parameters. Generally speaking, a ABIC of more than 2 (or 6) 
could be c onsidered as positive (or strong) evidence against a model. 

In [128'], based on the joint analysis of the SNIa+BAO+CMB data, Davis et al. systematically studied a number 
of dark energy models and made a comparison among them. Utilizing their AIG and BIG values, the "rank" of these 
models are determined. The AGDM model achieves the best of all the models due to its nice fit to the data and 
the economy of parameters. A series of models, including the constant w model, the GPL parametrization and the 
Gardessian expansion, can also provide comparably good fits but have more free parameters, thus less favored by the 
IGs. Moreover, the dynamical behavior of these models all tend to collapse to the AGDM. The DGP and the standard 
Ghaplygin gas models are clearly disfavore d. More studies on the comparison among dark energy models have been 
performed in with similar results obtained |l34l . Il36l . Il38l [l49l . [l60] . 

In all, although some dark energy models have been ruled out (or shown to be less optimistic) based on the current 
cosmological observations, we still need more observational data with higher precision to distinguish the various dark 
energy models. 

D. Model-independent dark energy reconstructions 

In addition to the cosmological constraints on specific dark energy models, another route of numerical studies, i.e. 
the model-independent dark energy reconstructions, has also drawn more and more attentions. 
The dark energy reconstruction is a classic statistical inverse problem for the Hubble parameter 

H{z) = HaV^nM^ + + (1 - ^mo)fiz), (31) 

where f{z) = Pde{z) / Pde{Q) is the dark energy density function. Different reconstruction method will give different 
/(z). The main target of a dark energy reconstruction is to detect the dynamical property of dark energy. It is widely 
believed that the EOS of dark energy Wde = Pde/ Pde holds essential clues for the nature of DE, because it is related 
with f{z) through an integration 

/(.) = exp(3jj^.'i±^). (32) 

The model-independent dark energy reconstructions can be divided into four classes: (i) Specific Ansatz: assuming 
a specific parameterized form for Wde{z) and estimating the associated parameters, (ii) Binned Parametrization: 
dividing the redshift range into different bins and using a simple local basis representation for Wde{z) or pde{z). (iii) 
Polynomial Fitting: treating the dark energy density function f{z) = Pde{z) / Pde{Q) as a free function of redshift and 
representing it by using the polynomial, (iv) Gaussian Process modeling: using a distribution over functions that can 
represent Wdeiz) and estimating the statistical properties thereof. These four classes of reconstruction methods and 
the related research works will be introduced in the following. 

1. Specific ansatz 

The "specific ansatz" is the most popular approach currently. The key idea is assuming a specific parameterized 
form for Wde{z) and estimating the associated parameters. A simple and widely used ansatz is the XGDM ansatz, in 
which the EOS of dark energy is a constant, i.e. Wde = const. This yields a simple form of f[z) 

/(z) = (l + z)3(i+"'«=). (33) 

In jl6l| . by combining the WMAP5 observations with BAO and SN data, Komatsu et al. obtained Wde = ~0.992lg'Qg2 
at the 1(7 GL, while in jl23| . a combination analysis of WMAP7-t-BA0-|-SN gave Wde = -0.980^^;!^^i^. A more recent 
constraint on Wde by the SGP team |88j| also presented the consistent result. So the current observations still favor 
Wde = — 1 (i-e. cosmological constant). 
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Besides, one can also assume that the EOS of dark energy is not a constant. The most popular parametrization 
with dynamical Wde is the Chevallier-Polarski-Linder (GPL) p arameterization jl62l . Il63| . which assume Wddz) = 
Wo + Waz/{1 + z). The corresponding /(z) is given by |l62l . Il64| 

/(z) = (1 + z)3(i+-o+-.) exp . (34) 

Here denotes the value of the present EOS, while Wa denotes t he va riation of the EOS. The current observational 
data, such as the Union2 SNIa dataset ^] and the WMAP7 data [123{ . still favor the result of = — 1 and Wa — 0, 
which is consistent with the cosmological constant. 

In addition, some other ansatzs have also been proposed, such as Wde{z) = wq+wiZ [165| |. Wde{z) = wo+Waz/{l+zy 
|l66j , wh^(z) = wn + wi ln(l + z) [l67j |. and Wdeiz) = wq + wi ( '"l^^-* — In 2) jl68j |. For more research works of various 
parametrization forms, see jl69l - ll72| and references therein. 



2. Binned parametrization 

In addition to the specific ansatz, another popular app roach is the binned parametrization. The binned parametriza- 
tion was firstly proposed by Huterer and Starkman |l73j | based on the principal component analysis (PGA) jl73l Il74| . 
The key idea is dividing the redshift range into different bins and picking a simple local basis repres enta tion for Wde{z). 
The simplest way is setting Wde{z) as piecewise constant in redshift, thus /(z) can be written as jl75l , [l76,] 

f{zn-i <z<z^) = {l + z)3(i+-") [](1 + z,f^^^-^'+^\ (35) 

i=0 

where Wi is the EOS parameter in the i-th redshift bin defi ned by an upper boundary at z^. 

The optimal choice of redshift bins is still in debate. In [l77| . Riess et al. proposed an uniform, unbiased binning 
method, in which the number of SNIa in each bin times the width of each bin is a constant (i.e. rtAz — const). In 
|l78j . Wang argued that one should choose a constant Az for redshift slices, to ensure the variation of H and \/Da 
in each redshift slice remain roughly constant with z. In [l79j |. we presented a new binned parametrization method, 
which treats Zi as models parameters and let them run freely in t he redshift region of SNIa samples. This binning 
method can obtain smaller Xmm and tighter error bars [l80l Il8l( . It should be mentioned that, since the current 
observational data at high redshifts (i.e.. z > 1) are very rare, for the binned parametrization, only two parameters 
of EOS Wde can be constrained well [l82j | . 

In addition to the piecewise constant parametrization, some oth er local basis representations for Wde{z) are also 
proposed, such as wavelet jl83l | and numerical derivatives |l84ll85| . 



3. Polynomial fitting 

The third approach is the polynomial fitting method. The key idea is treating the dark energy density function /(z) 
as a free function of redshift and representing it by using the polynomial. Gompared with the binned parametriza- 
tion, the advantage of the polynomial fitting parametrization is that the dark energy density function /(z) can be 
reconstructed as a continuous function in the redshift range covered by the observational data. 

A simple polynomial fit to /(z) was proposed by Alam et al. (186| , which is a truncated Taylor expansion 

f{z)^Ao + Ai{l + z)+A2{l + zf. (36) 
This ansatz has only three free parameters {^mO^ ^i, ^2) since Aq + Ay + A2 = 1 — ilmo for a flat Universe. Using this 



ansatz, Alam et al. 187l | argued that the Tonry/Barris SNIa sample jl88l | appear to favor dark energy whic h evo lves 
in time ^ . Another interesting polynomial fit is the polynomial interpolation, which was proposed by Wang |l9l| . It 
choose different redshift points Zi — i * z^i^x/n{i = 1, 2, • • • , n), and interpolate /(z) by using its own values at these 
redshift points. This yields 

^ (Zi ~ Zl)...[Zi - Zi^l){Zi - Zi+l)...[Zi - z„) 



It should be mentioned that there was a debate about the reUabiUty of this ansatz |189| . llQOH 
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Here fi = f{zi) and Zn = ^max- Based on the relation /(O) = 1, one parameter can be fixed directly, and only 
n — 1 model para meters need to be determined by the data. For more research works of this polynomial interpolation 
method, see Il92ll. 



4- Gaussian process modeling 

The fourth approach is the Gaussian Process (GP) modeling, which is proposed by Holsclaw et al. jl93l . llQ^ - GP 
is a stochastic process indexed by z. The defining property of a GP is that the vect or that corresponds to the process 
at any finite collection of points follows a multivariate Gaussian distribution fl95']. GPs are elements of an infinite 
dimensional space, and can be used as the basis for a nonparametric reconstruction method. T hey are characterized 
by a mean and a covariance function, defined by a small number of hyperparameters [l95l Il96j | . 

Based on the definition of a GP, one can assume that, for any collection zi, z„, Wdeizi), Wde(zn) follows a 
multivariate Gaussian distribution with a constant negative mean and exponential covariance function written as 

K{z,z') = t^x\''-''\\ (38) 

The hyperparameters x £ (0, 1) and k, and the parameters defining the likelihood, are determined by the data. 
The value of a G (0, 2] influences the smoothness of the GP realizations: for a = 2, the realizations are smooth 
with inflnitcly many derivatives, while a = 1 leads to rougher realizations suited to modeling continuous non-(mean- 
squared)-differentiable functions. Moreover, one can set up the following GP for Wde 

Wde{u)^GP{-l,K{u,u')). (39) 

Making use of the Eq. (P^, one ca n take advantage of the particular integral structure of luminosity distance di^iz) 
expressed by Wde{z) (see (193l . ll94| for details). 



IV. CONCLUDING REMARKS 



We have reviewed theoretical models as well as observational and numerical studies of dark energy. Numerous 
works and papers have been done and written since the discovery of the accelerating expansion of the Universe, it is 
impossible to cover even a small part of the heroic endeavors of our community in any review article. 

However, the problem of understanding the nature of dark energy is as daunting as ever, or perhaps some already 
hold the key to this understanding without being commonly accepted yet. Clearly, there is a long long way to go for 
both theorists and experimentalists. 

It is without any doubt that the process of detecting the nature of dark energy and understanding its origin will 
prove to be one of the most exciting stories in modern science. 
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